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FOREWORD 


This  work  was  carried  out  in  response  to  the  results  of  the  systems  analysis  conducted 
in  connection  with  the  preparation  of  the  QMDO  for  Functional  Field  Shelters.  This 
analysis  identified  the  pressurized  rib  concept  as  being  the  most  promising  structural 
alternative  with  regard  to  meeting  die  Army  requirement  for  lightweight  highly  mobile 
tentage.  A project  was  initiated  in  FY  71  to  investigate  the  feasibility  of  the  pressurized 
rib  concept  and  to  develop  the  design  data  on  pressure  stabilized  structural  elements 
necessary  for  Army  tentage  applications.  Support  for  this  project  has  been  provided 
through  the  In-House  Laboratory  Independent  Research  program. 

This  report  presents  a mathematical  analysis  of  the  deformation  behavior  of  these 
structural  elements  under  load  and  the  development  of  the  equations  which  will  be  used 
in  the  design  of  structures.  Subsequent  reports  will  present  experimental  studies  designed 
to  test  this  theoretical  analysis  and  will  describe  extensions  of  this  study  to  the  performance 
of  arches. 
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INTRODUCTION 


The  use  of  pressure  stabilized  structures  for  Army  field  shelters  has  been  the  subject 
of  considerable  attention  over  the  past  decade  with  emphasis  on  the  single-wall  and 
double-wall  type  shelters,  as  indicated  by  references  1 and  2.  A concept  that  has  received 
little  attention  is  that  of  constructing  tent  frames  with  pressure  stabilized  beams  and  arches. 
In  view  of  this,  an  investigation  of  the  behavior  of  these  structural  elements  under  load 
was  initiated  and  this  report  presents  the  portion  of  this  work  relating  to  the  deformation 
behavior  of  the  pressurized  beams.  Specifically,  a theoretical  study  of  the  problem  is 
presented  with  the  objective  being  the  establishment  of  a set  of  equations  to  be  used 
in  the  design  of  tent  frame  structures  using  pressurized  beams. 

The  structures  under  consideration  can  be  described  as  cylindrical  membranes,  and 
considerable  theoretical  work  has  been  done  on  such  structures.  The  linear  work  in  this 
area  is  a statically  determinate  problem  with  the  associated  limitations  regarding  the  types 
of  boundary  condition  that  can  be  applied.  These  limitations  are  too  restrictive  for  the 
problem  under  consideration  and  the  removal  of  these  restrictions  requires  the  use  of 
a nonlinear  membrane  theory.  Since  the  objective  of  this  work  was  development  of  a 
linear  theory,  the  classical  membrane  approach  was  rejected.  Most  other  work  regarding 
the  behavior  of  this  type  of  structure  can  roughly  be  divided  into  two  categories.  The 
first  of  these  included  work  which  examined  special  cases  such  as  beams  subjected  to 
uniform  moments  or  statically  determinate  loadings  in  which  the  stresses  could  be 
calculated  without  regard  to  deformation.  In  general,  the  work  in  this  first  category, 
e.g.  references  3 and  4,  does  not  constitute  attempts  to  develop  a generally  applicable 
theory  to  be  used  for  the  design  of  pressure  stabilized  beams.  The  second  category,  e.g. 
references  5,  6 and  7,  deals  with  the  wrinkling  behavior  of  pressurized  cylinders  which 
relates  to  the  failure  criteria  for  this  type  of  structure.  This  is  an  extremely  complicated 
problem  and  these  papers  therefore  deal  with  simplified  loading  problems,  typically  uniform 
moments.  Thus  these  results,  although  useful  in  developing  failure  criteria,  are  not 
applicable  to  general  loadings  and  boundary  conditions.  A theory  for  the  behavior  of 
pressure  stabilized  beams  is  presented  in  reference  8;  however,  this  theory  was  not  used 
in  the  present  work  because  it  was  developed  by  writing  expressions  for  the  3 components 
of  energy  independently  and  summing  them  to  obtain  the  total  energy.  It  was  felt  that 
a more  consistent  derivation  would  result  from  application  of  the  ideas  of  superposition 
of  small  displacements  on  large  ones,  that  is,  those  due  to  pressurization.  In  addition, 
the  results  in  reference  8 deal  principally  with  the  buckling  of  columns. 

The  present  report  includes  a derivation  of  the  governing  equations  and  boundary 
conditions  for  the  lateral  deformation  of  pressure  stabilized  beams  and  development  of 
Green's  Function  and  uniform  load  solutions  for  the  simply  supported  and  clamped  beam 
cases.  In  addition,  parametric  studies  are  presented  which  illustrate  the  linear  behavior 
of  pressure  stabilized  beams  under  load. 


ANALYSIS 


In  this  section  a derivation  of  a set  of  governing  equations  for  lateral  deformation 
of  pressure  stabilized  beams  and  their  solution  will  be  carried  out  The  derivation  will 
utilize  the  principle  of  minimum  potential  energy  and  begins  by  establishing  an  energy 
principle  for  deformation  about  the  pressurized  state.  A simplified  displacement 
approximation  is  then  used  to  reduce  the  problem  from  two  dimensions  to  one  dimension. 
This  displacement  approximation  requires  that  the  cross  section  remain  underformed  and 
includes  deformation  due  to  transverse  shear.  The  governing  equations  based  on  this 
approximation  are  then  obtained  using  variational  principles.  The  fundamental  homogenous 
solution  to  these  equations  is  obtained  in  preparation  for  development  of  the  Green's 
function  solutions  for  the  cases  of  simply-supported  ends  and  clamped  ends.  In  addition 
to  these  Green's  function  solutions,  the  solutions  for  the  uniformly  loaded  beams  are 
obtained  directly  by  finding  the  appropriate  particular  solution  of  the  governing  equations. 

Derivation  of  the  Governing  Equations 

Fundamental  Principles: 

The  derivation  is  based  on  the  principle  of  minimum  potential  energy  which  for  the 
membrane  state  of  stress  can  be  written: 

6/^  -iQiiN,!  + CjjNjj  + 2£]  j N i j — 2pv3J  dA  = 0 (1) 

The  subscripts  denote  the  coordinated  directions  Xj  and  x3  as  shown  in  Figure  1.  In 
equation  (1)  the  e;j  and  Njj  are  respectively  the  surface  strains  and  the  inplane  stress 
resultants;  p is  the  normal  load  distribution  and  v3  is  the  normal  displacement.  For  the 
displacements  illustrated  in  Figure  1 the  strains  are  defined  as 

en  = v,  , + 1 (v3,)2  (2a) 

e22  = — (v2  2 + v3)  + |-v3  2 + v2)2  (2b) 

r 2r2  ’ 


€12  * — [*2,1  + — V|  2 +—  (-V3(i)(-V32  + V2 )) 
2 r ' r 


where  the  commas  denote  differentiation  with  respect  to  the  coordinate  designated  by 
the  subscript  following  the  comma.  The  strei  resultants  are  given  in  terms  of  the  strains 
by  the  stress-strain  relation  as: 
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FIGURE  1.  COORDINATES  AND  DISPLACEMENT 
COMPONENTS. 
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(3c) 

Energy  Principle  for  Deformation  about  Pressurized  State: 

To  derive  the  energy  principle  for  deformation  about  the  pressurized  state  it  is  assumed 
that  the  displacement  components  Vj  can  be  written  as  a superposition  of  displacement 
components,  w,,  of  the  pressurized  state  and  Uj,  describing  the  deformation  about  the 
pressurized  state  due  to  applied  loads  as 

Vj  = Wj  + Uj  (4) 

Substituting  equation  (4)  into  the  strain  displacement  relation  (2)  yields 
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where 

Ci°i  * w,  i + 1 (w3>1)J  (6a) 

6i'i  = ulfl  + w3<1  u3fl  (6b) 

t‘i  i = -i-  (u3i)J  (6c) 

2 

6j°2  = — (Wj  j + W3 ) +_L  (— W3  j + Wj)J  (6d) 

r 2rJ 

£22  = J_  (Uj>3  + U3)  + JL  (-W3<2  + Wj)(-U3<j  + Uj)  (6e) 
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2r2 
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*12  = y [u2,,  + -1  U1<2  + 1 (-w3)1)(-u3|2  + u2) 

+ -1  {— U3  ,)(-W3  2 + W2)] 
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(6h) 


e'i2  = A-  [-u3(,(-u32  + u2)l 

By  the  use  of  equations  (5)  in  the  stress-strain  law  (3)  obtain 
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Substitution  of  equations  (5)  and  (7)  into  the  energy  principle  and  retaining  terms  up 
to  those  quadratic  in  the  displacement  components  u;  yields 

6 / a _1~ l(fii  N , , + e2 2 Nji  + 2eia  Nj  2 — 2pw3 ) + 2 (ei  j N j j 

« 2 

+ ej j Nj j + 2eu  N, j — PU3)  + (cii  N, , + ej 2 Nj j + 2e, 2 Nj2 


+ 2 £ 11  M11  + 2e2  2 N22  + 4ei2  N12  — 2FU3)]  dA  = 0 


(9) 


In  establishing  this  relationship  the  energy  due  ♦ o the  applied  forces  F which  do  work 
only  through  the  displacement  component  U3  has  been  added  and  the  following  identities 
which  can  be  proven  through  the  stress-strain  law  have  been  used: 
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To  yield  the  governing  equations  for  deformation  aoout  the  pressurized  state,  which 
is  described  by  the  displacement  components  Uj,  the  variation  in  equation  (9)  must  be 
interpreted  so  as  to  be  carried  out  with  respect  to  the  components  Uj.  Since  the  expression 
inclosed  in  the  first  set  of  brackets  in  equation  (9)  is  independent  of  the  components 
u;,  its  variation  must  vanish.  Carrying  out  the  variation  of  the  second  set  of  brackets 
and  performing  the  required  integration  by  parts,  one  finds  that  the  coefficient  of  the 
variational  displacements,  5u:,  are  the  equilibrium  equations  governing  tSie  pressurized  state 

' 0 0 o 

which  must  vanish  if  the  stress  resultants  N,,.  N22  and  N,2  represent  an  equilibrium 
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state.  Taking  this  into  account,  the  following  is  the  energy  principle  for  deformation 
about  the  pressurized  state: 

8 f ^ — (cj t N , | + ejj  N12  + 2en  Njj  + 2ei  i Nj  i + 2ejj  N21 

2 

+ 4eJ2  N,j  - 2FUj)  dA  = 0 (11) 

This  is  the  energy  principle  to  be  used  to  derive  the  equations  governing  the  deformation 
about  the  pressurized  state. 

Simplified  Displacement  Approximation: 

The  problem  as  presently  stated  in  equations  (11),  (6)  and  (8)  is  a problem  in  two 
independent  coordinates  and  will  yield  partial  differential  equations.  The  objective  here 
if  to  simplify  the  mathematical  form  by  introducing  a displacement  approximation  which 
specified  the  form  of  the  displacement  components  with  respect  to  the  circumferential 
coordinate.  This  displacement  approximation  describes  the  deformation  in  terms  of  rigid 
body  motions  of  the  cross  section  or  stated  alternatively,  constrains  the  cross  section  to 
retain  its  shape  during  deformation.  This  approximation  is  the  one  used  in  beam  theory; 
here,  however,  transverse  shear  deformation  is  included.  This  approximation  can  be  stated 
mathematically  by  expressing  the  displacement  components  U;  in  terms  of  the  axial 
displacement  of  the  cross  section,  U,  the  tranwerse  displacement  of  the  cross  section, 
W,  and  the  rotation  of  the  cross  section  about  an  axis  normal  to  the  plane  of  0 and 

73,  <t>,  as: 

U]  = U - <t>a  Sin(x2)  (12a) 

u2  =WCos(Xi)  (12b) 

uj  = WSin(xj)  (12c) 


These  cross  section  displacements  which  are  functions  of  x 1 only  are  illustrated  in  Figure  1. 
These  displacement  approximations  can  be  used  to  formulate  one-dimensional  strain 
measures,  stress  resultants  and  a one-dimensional  energy  principle.  The  one-dimensional 
strain  measures  are  obtained  by  substituting  equations  (12)  into  equations  (6).  This  yields 
the  following  non-vanishing  strain  component;: 
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ej , * e,  - x,  a Sin(xj)  (13a) 

c”,  «*  J_  fiR’  )2  Sin(x2 ) (13b) 

2 

c,’2  « e,2  Cos(x2)  (13c) 

In  carrying  out  this  formulation  the  displacement  components  due  to  pressurization  w2 
and  w3  have  been  assumed  to  vanish  and  take  on  a constant  value,  respectively.  The 
one-dimensional  strain  measures  e1(  ku  and  e13  are  defined  as: 

e,  = 1 U'  (14a) 

P 

K,  = -!*'  (14b) 

6,2  = 1 (1  W ' - 0)  (14c) 

2 p 

In  writing  these  expressions  the  axial  coordinate  and  the  displacements  have  been 
non-dimensionalized  as  follows: 

x>  - l{  (15a) 

W * aW  (15b) 

U « aU  (15c) 

and  the  prime  denotes  differentiation  with  respect  to  the  non-dimensional  variable,  $,  and 

a and  I are  respectively  the  beam  cross  sectional  radius  and  length  measure. 

The  one-dimension»!  energy  principle  is  established  by  substitution  of  equations  (13) 
into  equation  (11),  taking  account  of  the  vanishing  strains  this  yields: 

5 1 [e,r,,  - K,n( , aSin(x3)  + 2eiania  Cos(x3) 

+ _L(W'iJn  Sin2(x3)  -2  1 F W Sin(x2)  I aldx3df  - 0 (16) 

P1  P 
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In  this  expression  the  stress  resultants  have  been  non-dimensionalized  by  division  by  a 
reference  value  of  the  elastic  modulus  Cn  = C, The  non-dimensional  force  parameter 
is  defined  by  F = Fl/C, , . This  energy  principle  is  reduced  to  a one-dimensional  principle 
by  definition  1/  the  o';e-dimensional  stress  resultants  and  forces  as  follows: 


T*//  nn  adxi 

(17a) 

M=  n, , a2Sin(x2)dx2 

(17b) 

2 7T 

Q = f0  n,2  aCos(x2)dx2 

(17c) 

f = - sV  FSin(x2)dx2 

(17d) 

The  remaining  term  appearing  in  the  energy,  that  containing  the  stress  due  to  pressurization, 
n,  can  be  integrated  with  respect  to  x2  since  both  n and  W are  functions  cf  x,  only. 
The  one-dimensional  stress-strain  law  is  obtained  from  equations  (17)  by  substitution  of 
the  stress  strain  law  (8),  properly  non-dimensionalized,  and  the  strain  displacement 
relations  (13)  and  carrying  out  the  required  integration  to  yield 


T * 27rdae, 

M = -irda3  k, 
Q = 2?rcae,  7 


(18a) 

(18b) 

(18c) 


The  one-dimensional  energy  principle  is  put  in  final  form  by  rewriting  equation  (16)  using 
tne  definitions  given  in  (17): 


6 / — [e,T  - k,M  + 2ej 2 Q + **  (W')2n 

2 Pr 


- 2tt  1 fW]  Id*  = 0 
P 


(19/ 


It  should  be  remembered  that  the  sheas  resultant  due  to  pressurization,  n,  is  a known 
quantity  in  this  equation. 
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Governing  Differential  Equations: 


With  the  preliminaries  as  developed  above,  it  is  possible  to  derive  the  equations 
governing  the  cross  section  displacements  U,  W and  4>-  This  is  accomplished  by  writing 
the  energy  principle  (19)  in  terms  of  these  displacements  by  using  equations  (18)  and 
(14)  and  carrying  out  tiie  variation,  die  integration  by  parts  where  required  and  application 
of  the  fundamental  lemma  of  the  calculus  of  variation  to  yield  the  differential  equations. 
The  energy  principle  in  terms  of  displacements  is  given  as: 


* ft  -^L  (U')J  + 4-  l*')2  + c(-  W-6)1  + -1  (W')J 


2 pJ 


-2 — fW)  d£  = 0 (20) 

P 


carrying  out  the  variation 


ft  ?ral  (-4  U'5U'  + ~ $'6  0'  + W'6  W*-  1 W'5£  - - <t> 6W' 

* P p 2 pi  p p 

+ c060  + JL  W'6W*  - 1 f6W)  d$  = 0 (21) 

P2  P 

Carrying  out  the  required  integration  by  parts  with  the  interval  of  integration  being  (£,, 

i 2 

/ na\  ((_  id  u-.,  8U  + (_  d _ c w'  + o^)  54, 

*1  Di  Dl  n 


+ (-  1 (C  + n)W”  + £ 0'  - 1 f)  6W]  d{ 
P2  p p 


+ [1  (c  + n)W'  -i  + (2d  y.j  6U/fc  + jd  60/fc 

A*  A f • A*  >•  A > . 


1 p* 


IJ  p 


According  to  the  fundamental  lemma  of  the  calculus  of  variations  satisfaction  of  this 
equation  requires  that  the  coefficients  of  £U,  (W  and  vanish  independently  yielding 
the  governing  differential  equations 


U"  = 0 


(23) 


d0"  + cpW'  - cp*0  = 0 


(c  + n)W"  - cp0'  + pf  = 0 


and  boundary  conditions  requiring  specification  of 


U‘  or  U 


0 or  0 


1 (c  + n)W'  - i 0 or  W 
P2  P 


at  Mi  and  Ma- 
lt is  to  be  noted  that  in  these  euqations  the  axial  deformation  specified  by  the 
dependent  variable  U is  uncoupled  from  the  trc.<sverse  deformation  specified  by  the 
dependent  variables  W and  0.  Thus  when  axial  deformation  is  of  no  interest,  as  is  the 
present  case,  it  may  be  ignored. 

Solution  of  the  Governing  Equations 

The  purpose  of  this  section  is  the  formulation  of  the  solutions  to  the  equation 
obtained  in  the  previous  section.  This  work  will  be  limited  to  the  transverse  or  bending 
deformation  parameters,  W and  0.  The  fundamental  homogenous  solutions  will  be  obtained 
first  and  this  will  be  utilized  to  develop  Green's  function  solutions  for  the  cases  of  simply 
supported  ends  and  clamped  ends. 

Fundamental  Homogenous  Solution: 

In  this  section  the  homogenous  solution  to  the  equations  governing  W and  0 are  found. 
These  equations,  (24)  and  (25),  are  second  order  coupled  ordinary  differential  equations 
with  constant  coefficients  so  die  solution  is  of  the  form: 


I 


W * Aew£ 


(27a) 


0 = Bew$ 


Substitution  of  these  expressions  into  the  homogenous  differential  equations  leads  to  the 
following  matrix  equations  for  the  constants  A and  B and  the  characteristic  roots 


(c  + n)to2  -pccJlA 


dor  -cp  B 


For  these  equations  to  have  a nontrivial  solution  the  determinate  must  vanish  from  which 
the  characteristic  roots  are  found  to  be  u>i  = o,«j  = 0,  u3  = X,  and  u3  = -X.  For 
each  of  these  roots  a relation  exists  between  A and  B,  as  follows: 

e = A ,/p 
B,  a 0 

03  - 0*4  (29) 


04  * 0A3 


p = cpV(cp2  - X2d) 


Vd(c  + n) 


Thus  the  homogenous  solution  is 


W * A,  +A2$  + A3  Cosh(X£)  + A4  Sinh(X^) 


0 = •*—  Aj  + A 40  Cosh(X{-)  + A3P  Sinh{X$) 
P 


Green's  Function: 


The  solution  obtained  in  the  previous  section  is  not  complete  in  that  neither  the 
loading  condition  nor  the  boundary  conditions  have  been  specified.  In  this  section  two 
sets  of  boundary  conditions  will  be  considered  in  conjunction  with  a concentrated  loading. 
This  will  be  handled  in  the  form  of  a Green's  function  because  the  results  are  applicable 
not  only  to  the  concentrated  load  but  also  to  the  general  loading  situation  by  means 
of  direct  integration.  To  carry  out  the  solution,  consider  a beam  of  length  21  with 
the  origin  of  coordinates  located  at  the  center  of  the  beam  and  concentrated  load,  g, 
applied  at  $ = tj.  In  a physical  situation  this  load  would  be  applied  as  a distributed 
circumferential  line  load  G(x2)  and  this  is  related  to  the  concentrated  load,  g,  by: 

1 21t 

g = — /0  G(x2)  Sin(x2)dx2  (31) 

irC, , 

The  Green's  function  is  found  (reference  9)  by  obtaining  two  sets  of  functions  of  the 
form  of  equations  (30),  one  of  which  (W, , 0, ),  satisfies  the  boundary  conditions  at  $ = -1 
and  the  other,  (W2,  02),  satisfies  the  boundary  conditions  at  £ * 1.  The  remaining 
constants  of  integration  are  determined  by  satisfaction  of  three  continuity  conditions  and 
a jump  condition  to  account  for  the  concentrated  load  at  £ * rj.  These  conditions  are: 

W.fo)  « W,to)  (32a) 

0i  (t?)  = (t?)  (32b) 

0i(f?)  = 02  (t?)  (32c) 

[(c  + n)W j(r?)  - cp0,(r?)j  - [(c  + n)W2'(?j)  - cp02(t?)l  = P9  . (32d) 

Simply  Supported  Boundary  Conditions: 

The  first  case  of  end  restraint  to  be  considered  is  that  of  a beam  simply  supported 
on  both  ends.  Physically,  simply  supported  means  that  the  transverse  displacement  and 
the  moment  vanish  at  the  end  of  the  beam.  These  conditions  are  expressed  mathematically 
as 
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w,  (-1)  - 0 


0l't— 1)  = o 


Wa(1)  = 0 


0,'(1)  » 0 


Application  of  these  four  conditions  in  conjunction  with  the  four  conditions  at  the  location 
of  the  concentrated  load,  equations  (32),  is  sufficient  to  determine  the  eight  constants 
appearing  in  W1(  0,,  W2  and  02,  thus  completing  the  solution  of  this  case  which  is  as 
follows: 


for  -1  < E < t? 


W = JL  (p(1-tj)(1+£)  + r,  (Sinh(XE)  + Tanh(X)  Cosh(Xf)l) 
2n 


0 = -JL  [(1-tj)  + /3r,  [Cosh(XE  )<■  Tanh(X)  Sinh(Xf))) 


T * Sinh(Xq)  - Tanh(X)  CosMXt?) 
0Tanh(X) 


for  t?  < { < 1 


W «JL  (p(1+t})(1-{)  + r2  [Sinh(XE)  - Tanh(X)  Cosh(Xf)]] 
2n 


0 = i L l— (1—7?)  f 0 r2  (Cosh(XE)  - Tanh(X)  Sinh(XE))] 

2n 


P _ Sinh(Xq)  + Tanh(X)  Cosh(XT?) 
/JTanh(X) 


In  writing  these  expressions  the  subscripts  on  W and  0 have  been  dropped  and  the  range 
of  E over  which  the  solution  is  valid  has  been  specified.  This  solution  has  the  symmetry 
property  required  of  a Green's  function,  namely  W(E,r?)  - W(»?,E). 


■ 


In  this  section  the  Green's  function  for  a beam  with  clamped  or  fixed  restraints 


on  both  ends  will  be  determined.  Physically,  the  clamped  end  means  that  both  the 
transverse  displacement  and  the  slope  vanish  at  the  end  of  the  beam.  This  is  expressed 
mathematically  as 

W,  (-1)  = 0 

(38a) 

o 

ii 
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(38b) 

Wj(1)  = 0 

(38c) 

0j( 1)  = 0 

(38d) 

As  with  the  simple  supported  case,  application  of  these  conditions  in  conjunction  with 
t those  specifying  the  condition  at  the  position  of  the  concentrated  load  allows  all  unknown 

constants  of  integration  to  be  determined  giving  the  solution  as 

for  -1  > £ > r/ 

W *1_  [a3[p/J(1+£)  Sinh(X)  - Cosh(X)  + Cosh(X£)] 

2n0  (39a) 

> + a4[-p/3(1+£)  Cosh(X)  + Sinh(X)  + Sinh(X£)]] 

0 « JL  [a3  (Sin(X)  + Sinh(X£)]  + a4  (-Cosh(X)  + Cosh(X£))  (39b) 

2n 

for  17  > £ > 1 

W =2_  [b3  [p0(1— £)  Sinh(X)  -Cosh(X)  + Cosh(X£)J 
2nj3 

+ b4  [p0(1— £)  Cosh(X)  — Sinh(X)  + Sinh(X£))J  (40a) 

[b3[— Sinh(X)  + Sinh(X£)]  + b4  (-Cosh(X)  + Cosh(X£)]  ] (40b) 

in 

15 


I 


where 


a3  * (1  + Sinh(X)  SinhIXr?)  - Cosh<X)  Cosh(Xr})]  /Sinh(X)  (41a) 

84  * [-pPn  - pP  Sinh(X)  Sinh(Xij)  + p0  Cosh(X)  Cosh(Xr?) 

(41b) 

+ Cosh(X)  Sinh(Xtj)  - Sinh(X)  Cosh(Xr?))/(Sinh(X)  -pp  Cosh(X)) 
b3  = (1-  Sinh(X)  Sinh(XTj)  - CosH(a,  Cosh(r?X)] /Sinh(X)  (42a) 

b4  = 1~pPv  - PP  Sinh(X)  Sinh(Xr?)  - p/?Cosh(X)  Cosh(Xr)) 

(42b) 

+ Cosh(X)  Sinh(Xr?)  + Sinh(X)  Cosh(Xr?))]  /|Sinh(X)  - pp  Cosh(X)) 


t 

I 


L 


Although  it  cannot  be  seen  by  inspection,  this  solution  also  has  the  symmetry  property 
associated  with  a Green's  function. 


Uniform  Load  Solution: 

A solution  that  is  frequently  useful  and  relatively  easy  to  obtain  directly  as  opposed 
to  integration  of  the  Green's  function  is  that  for  the  uniformly  loaded  beam.  As  with 
the  Green's  function,  solutions  will  consider  a beam  of  length  21  with  the  origin  of 
coordinates  at  the  center  of  the  beam  as  depicted  in  Figure  V The  solution  desired 
is  that  for  equations  (24)  and  (25)  with  the  loading  parameter,  f,  taken  as  a constant. 
This  solution  can  be  written  as  the  sum  of  a homogenous  solution  and  a particular  solution. 
The  homogenous  solution  is  given  by  equations  (30)  so  only  the  particular  solution  remains 
to  be  determined  and  this  is  taken  in  the  form: 

W = D|  (-1  (43a) 

<t>  • D2*  (43b) 

Substitution  of  these  into  the  differential  equations  (24)  and  (25)  and  solution  for  the 
constants  D,  and  Dj  yields 

D,  = -Jl  (44a) 

2n 
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D,  = - 1 (44b) 

n 

and  the  full  solution  is  then: 

W = I [-  h1  + A,  + As£  + A,  Cosh(X£>  f A4  Sinh(X£)]  (45a) 

n 2 

<t>  = l[-£  + !a,  + A40  Cosh(X£)  + A,0  Sinh(Xf)]  (45b) 

n P 

The  constants  Aj  are  determined  by  specification  of  the  boundary  conditions  at  £ = 1 
and  £ = -1. 

Simply  Supported  Boundary  Conditions: 


Application  of  conditions  (33)  in  conjunction  with  this  solution  (45)  gives  the  solution 
for  the  uniformly  loaded  beam  as  the  simply  supported  boundary  as 

W =1  [p(1+£J)/2  + [ Cosh (X£ ) — Co$h(X))] /X/JCosh(X)J  (46a) 

\ n 

4>  = i(-£  + Sinh(X£)ACosh(X)]  (46b) 

n 

, Clamped  Boundary  Conditions: 

Application  of  conditions  (38)  in  conjunction  with  the  solution  given  by  (45)  gives 
the  solution  for  the  uniformly  loaded  beam  with  clamped  boundary  conditions  as: 

W *1  [p(1— £2)/2  + (Cosh(X£)  - Cosh(X)]/0Sinh(X)]  (47a) 

n 

<p  - ! [-£  + Sinh(X£)/Sinh(X)l  (47b) 

n 

Two  observations  can  be  made:  1)  the  solutions  obtained  for  the  uniformly  loaded 

beam,  equations  (46)  and  (47)  are  identical  with  the  solution  that  would  be  obtained 
by  integrating  the  respective  Green's  functions;  2)  the  general  solutions  obtained  here 
(30)  and  (45)  can  be  used  to  obtain  solutions  for  other  combinations  of  boundary 
conditions  as  allowed  by  equations  (26). 
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DISCUSSION  OF  RESULTS 


The  purpose  of  this  section  is  the  presentation  of  results  which  describe  the  behavior 
of  pressure  stabilized  beams  under  load  as  predicted  by  the  theory  developed  in  the 
preceding  sections.  The  results  to  be  presented  will  describe  the  deformation  behavior 
and  the  strength  behavior  for  beams  subjected  to  a centrally  located  concentrated  load. 
This  problem  is  considered  here  because  it  represents  a characteristic  solution  of  the 
governing  equations  and  is  thus  very  useful  as  a tool  to  illustrate  the  behavior  with  respect 
to  the  three  parameters,  pressure,  geometry  and  material,  which  occur  in  the  equations. 
These  results  will  be  presented  for  two  cases  of  end  restraint,  simply  supported  and  fixed 
ends. 

Deformation  Behavior 

For  linear  structures  the  deformation  behavior  is  conveniently  characterized  by  the 
Green's  function  G(£,  17)  which  can  be  thought  of  as  the  deformation  at  a coordinate 
£ due  to  a unit  load  at  coordinate  17  or  the  flexibility  function  often  used  in  structural 
mechanics.  The  flexibility  parameter  used  here  is  the  deformation  at  the  beam  midpoir.t 
due  to  a unit  load  applied  at  the  same  point.  This  flexibility  parameter  is  presented 
in  nondimensional  form  as  a function  of  the  nondimensional  pressure,  geometry  and 
material  parameters. 

Pressure  Parameter: 

The  nondimensional  flexibility  parameter  as  a function  of  the  pressure  parameter 
is  shown  in  Figures  2 and  3 for  the  simply  supported  and  fixed  end  restraints,  respectively. 
For  the  data  shown  the  geometry  parameter  p = 100  and  the  modulus  ratio  is  c = 0.1. 
The  physical  pressure  range  covered  in  these  plots  depends  on  the  geometry  and  elastic 
modulus  of  the  particular  beam,  but  for  a beam  having  a radius  of  3.175  cm  (1.25  in) 
and  an  elastic  modulus  of  2.1015  x 10®  dyne/cm  (1200  Ib/in)  the  pressure  range  is 
0 - 4.236  x 106  dyne/cm2  (24.2  Ib/in2).  Examination  of  the  plots  in  Figures  2 and 
3 reveals  that  for  very  low  pressures  the  flexibility  is  large  and  rather  large  decreases 
in  flexibility  result  from  small  increases  in  pressure.  In  the  higher  pressure  range  the 
opposite  effect  is  observed,  that  is,  large  ircreases  in  pressure  are  required  to  significantly 
decrease  the  flexibility.  Thus,  in  the  higher  pressure  range  the  benefits  to  be  gained  with 
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ENSIONAL  FLEXIBILITY 


SION AL  FLEXIBILITY 


respect  to  decreasing  f lexit.  ility  or  increasing  stiffness  by  changes  in  pressures  are  not 
significant  Although  it  is  not  exactly  the  case  the  flexibility  can  generally  be  described 
as  varying  inversely  with  the  pressure.  The  same  type  of  behavior  is  observed  for  other 
geometries  and  modulus  ratios.  Compariso".  of  the  results  for  the  two  types  of  end 
restraints  reveals  that,  as  would  be  expected,  the  simply  supported  case  is  more  flexible 
than  the  fixed  case. 

Geometry  Parameters: 

Figures  4 and  5 show  respectively  the  variation  of  the  flexibility  parameter  with 
geometry  for  the  simply  supported  and  fixed-end  restraints.  As  indicated  above,  the 
geometry  is  specified  by  the  ratio  of  the  beam  length  to  the  cross-section  radius.  The 
results  presented  are  for  a modulus  ratio  c - G.1  and  for  three  values  of  the  pressure 
parameter.  These  results  reveal  a linear  relationship  between  flexibility  and  geometry 
except  in  the  low  range  of  the  geometry  parameter. 

Modulus  Ratio: 

The  third  parameter  appearing  in  the  governing  •’quations  is  the  ratio  of  the  shear 
modulus  to  the  elastic  modulus,  referred  to  hereafter  as  the  modulus  ratio.  The  flexibility 
parameter  as  a function  of  the  modulus  ratio  is  presented  in  Figures  6 and  7 for  the 
fixed-end  beam  and  in  Figures  8 and  9 for  the  simply  supported  beam.  For  each  case 
results  are  given  for  two  values  of  the  geometry  parameter  8nd  three  values  of  the  pressure 
parameter.  Examination  of  the  results  for  the  fixed-end  case  (Figures  6 and  7)  reveals 
that  for  short  beams  with  large  cross  sections  (Figure  6)  the  variation  of  the  modulus 
ratio  in  the  lower  range  of  that  parameter  can  have  a significant  effect  on  the  flexibility 
parameter.  This  effect  is  more  pronounced  at  low  pressures.  In  the  higher  range  of 
the  modulus  ratio  for  short  large-cross-section  beams  and  over  the  entire  range  of  modulus 
ratio  for  long  thin  beams,  the  flexibility  parameter  is  nearly  constant.  Thus,  it  is  apparent 
that  shear  deformation  is  an  important  consideration  for  short  beams  with  large  cross 
sections  utilizing  low-inflation  pressures.  Similar  results  are  observed  for  the  simply 
supported  case  (Figures  8 and  9)  although  the  importance  of  shear  deformation  is  not 
as  significant  as  in  the  fixed-end  case. 
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FIGURE  4.  FLEXIBILITY  AS  A FUNCTION  OF 

GEOMETRY,  SIMPLY  SUPPORTED  ENDS. 


22 


I 


NON  DIMENSIONAL  FLEXIBILITY 


0 01  0.2  0.3  0.4  03  0.6  0.7 

MODULUS  RATIO,  C 


FIGURE  6.  FLEXIBILITY  AS  A FUNCTION  OF 
MODULUS  RATIO,  FIXED  ENDS. 
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MODULUS  RATIO,  C 


FIGURE  7. 


FLEXIBILITY  AS  A FUNCTION  OF 
MODULUS  RATIO,  FIXED  ENDS. 


NONDIMENSIG 
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FIGURE  8.  FLEXIBILITY  AS  A FUNCTION  OF 
MODULUS  RATIO,  SIMPLY 
SUPPORTED  ENDS. 
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ON  DIMENSIONAL  FLEXIBILITY 


Failure  Criteria 


\ The  main  thrust  of  this  work  is  the  development  of  design  equations  for  use  on 

structural  elements  fabricated  of  materials  such  as  fabric  which  cannot  support  compressive 
stress.  For  this  situation  a useful  and  convenient  failure  criterion  for  describing  the  strength 
of  the  structural  element  is  the  wrinkling  load.  Several  authors,  see  for  example  references 
6 and  7,  have  carried  out  nonlinear  analysis  of  special  loading  situations  which  relate 
. this  wrinkling  load  to  the  ultimate  strength  of  the  member.  Thus,  this  wrinkling  load 

can  be  used  to  develop  other  failure  criteria.  This  wrinkling  load  is  defined  as  the  load 
which  results  in  a negative  stress  at  some  point  in  the  structure  equal  in  magnitude  to 
the  positive  stress  at  that  point  caused  by  pressure.  This  results  in  a stress  of  2ero  magnitude 
at  the  point  in  question  with  any  further  increase  in  the  load  causing  a compressive  stress 
which  cannot  be  supported  by  the  material  so  that  wrinkling  results.  The  wrinkling  load 
thus  specifies  the  initiation  of  failure  and  can  be  thought  of  as  analogous  to  the  yield 
j strength. 

As  was  done  with  the  discussion  of  the  deformation  behavior,  the  principal  results 
presented  here  will  be  for  the  case  of  a concentrated  load  applied  at  the  center  of  the 
beam  for  which  the  maximum  stress  and  thus  the  wrinkling  occurs  at  the  center  of  the 
beam.  T.esults  relative  to  the  behavior  of  the  wrinkling  load  as  a function  of  the  pressure, 
geometry  and  material  properties  will  be  presented  for  both  the  simply  supported  and 
\ fixed-boundary  restraints.  It  is  possible  to  express  in  a rather  concise  and  useful  form 

the  behavior  of  the  wrinkling  load  for  the  uniformly  loaded  beam  as  a function  of  the 
geometry  and  pressure  level  and  this  result  will  be  presented. 

Pressure  Parameter: 

» 

The  wrinkling  load  as  a function  of  pressure  level  is  depicted  in  Figure  10  for  the 
geometry  parameter  p = 100  and  the  modulus  ration  c * 0.1.  The  data  shown  is  for 
the  case  of  clamped  restraints  as  indicated,  but  the  data  for  the  case  of  simply  supported 
restraints  is  nearly  indistinguishable  from  the  results  presented  in  Figure  10.  Examination 
of  these  results  shows  that  the  wrinkling  load,  and  thus  the  load  supporting  capacity  of 
pressure  stabilized  beams,  increase  with  pressure.  The  rate  of  increase  is  greater  than 
a linear  relationship;  that  is,  doubling  the  pressure  results  in  more  than  a doubling  of 
the  wrinkling  load. 
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Geometry  Parameter: 


Figure  11  shows  the  behavior  of  the  wrinkling  load  as  a function  of  the  geometry 
parameter  for  both  the  clamped  and  simply  supported  end  restraints.  Data  is  shown  for 
three  values  of  the  pressure  parameter  n,  and  all  data  is  for  one  value  of  the  modulus 
ratio,  c = 0.1.  For  values  of  the  geometry  parameter  above  50  for  the  fixed  case  and 
30  for  the  simply  supported  case,  the  wrinkling  load  takes  on  a nearly  constant  value. 
Below  these  values  the  geometry  parameter  causes  considerable  variation  in  the  wrinkling 
load.  The  overall  variation  shown  in  Figure  11  is  that  of  the  negative  hyperbolic  tangent. 
These  resuits  show  that  above  some  value  of  the  geometry  parameter  the  load  carrying 
capacity  of  the  beam  is  independent  of  length.  Subsequently,  similar  behavior  will  be 
;hown  for  the  uniformly  loaded  beam. 

Modulus  Ratio: 

The  variation  of  the  wrinkling  load  with  the  modulus  ratio  is  shown  in  Figure  12 
for  the  case  of  simply  supported  end  restraints,  and  Figure  13  for  fixed-end  restraints. 
Results  are  shown  for  the  pressure  parameter,  n = 0.015  and  for  two  values  of  the 
geometry  parameter,  p = 20,  100.  These  results  show  that  for  values  of  c above  0.2 
the  modulus  ratio  has  very  iittle  effect  on  the  wrinkling  load  for  the  values  of  p and 
n presented.  These  results  are  typical  for  all  values  of  p and  n.  This  confirms  the  finding 
revealed  by  examination  of  the  deformation  behavior  with  regard  to  the  significance  of 
shear  deformation. 


Behavior  Under  Uniform  Load 


Since  the  uniformly  loaded  beam  more  closely  represents  the  type  of  loading  situation 
that  would  occur  in  tentage  application,  some  results  concerning  the  magnitude  of  the 
wrinkling  load  for  this  case  are  presented.  The  wrinkling  load  here  is  the  magnitude 
of  the  load  per  unit  length  of  beam  which  causes  wrinkling  to  begin.  The  results  presented 
here  are  for  the  load  per  unit  physical  length  as  opposed  to  the  load  per  unit 
nondimensional  length  used  in  the  derivation. 


For  simply  supported  end  restraints  the  wrinkling  load  is  given  as 


f 


w 


n2  Cosh  X 

a - Cosh  X + 1 


(48) 
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FIGURE  13 


WRINKLING  LOAD  AS  A FUNCTION 
OF  MODULUS  RATIO,  FIXED  ENDS 


The  parameter  X increases  with  both  n,  the  rji  assure  parameter  and  p the  geometry 
parameter  and  the  expressions  in  brackets  in  . nations  (48)  and  (49)  rapidly  approach 
unity  as  n and  p increase.  Thus,  over  much  of  the  pressure  and  geometry  range,  the 
wrinkling  load  for  both  types  of  end  rertraint  can  be  approximated  by 

fw  = £L  (50) 

2Cn 

when  a is  the  radius  of  the  tube  cross  section.  From  this  expression  it  is  seen  that  the 
load-carrying  capacity  increases  quadratically  with  the  pressure  and  linearly  with  the  cross 
section  radius.  Thus,  if  the  pressure  is  doubled  the  cross  section  radius  may  decrease 
by  a f'otor  of  four  and  the  same  load-carrying  capacity  is  maintained.  This  is  one  of 
the  principal  advantages  of  the  pressurized  rib  concept;  that  is,  the  use  of  highe-  pressure 

k 

levels  allows  a significant  decrease  in  the  cross  section  size. 

For  the  case  of  the  concentrated  load,  it  was  shown  that  above  some  values  of  the 
geometry  parametei  the  wrinkling  load  was  independent  of  the  geometry.  Examination 
of  equation  (50)  revoals  that  the  wrinkling  load  is  independent  of  the  beam  length  as 
► long  as  the  approximation  used  in  setting  the  terms  in  equations  (48)  and  (49)  to  unity 

are  valid. 

CONCLUDING  REMARKS 

f The  development  of  a linear  theory  for  the  behavior  of  pressure  stabilized  beams 

under  static  load  has  been  presented  and  solutions  to  the  resulting  governing  equations 
obtained.  Results,  principally  tor  the  case  of  concentrated  loading,  are  presented  which 
* illustrate  the  behavior  of  pressure  stabilized  beams  under  load  and  the  manner  in  which 

that  behavior  is  influenced  by  pressure  level,  beam  geometry  and  material  properties. 

With  regard  to  deformation  behavior,  the  results  show  that  the  flexibility  of  pressure 
stabilized  beains  has  a variation  with  pressure  that  is  very  nearly  an  inverse  relationship 
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and  a linear  variation  with  the  length  to  radius  ratio  except  for  small  values  of  that  ratio. 
The  material  properties  enter  in  two  ways,  first  through  the  pressure  parameter  which 
is  the  ratio  of  the  axial  stress  due  to  pressure  to  the  elastic  modulus,  and  sec*  nd  through 
the  modulus  ratio,  the  ratio  of  the  shear  modulus  to  the  elastic  modulus.  Thus,  the 
elastic  modulus  has  a direct  bearing  on  the  behavior  by  the  manner  in  which  >t  determines 
the  pressure  parameter.  The  results  show  that  the  influence  of  the  shear  modulus  on 
the  deformation  behavior  is  significant  only  when  modulus  ratio  is  less  than  0.2  and  then 
only  for  low  pressure  levels  and  length  to  radius  ratio. 

The  results  show  that  the  wrinkling  load  which  is  used  to  characterize  the  load-carrying 
capacity  of  pressure  stabilized  beams  increases  with  pressure  at  a rate  more  rapidly  than 
a linear  <elationship.  The  relationship  between  the  wrinkling  load  and  the  length  to  radius 
ratio  useu  to  (Ascribe  the  geometry  is  that  of  a negative  hyperbolic  tangent.  As  a 
consequence,  the  wrinkling  load  is  independent  of  the  geometry  for  values  of  the  length 
to  radius  ratio  greater  than  50.  As  with  the  flexibility,  it  is  found  that  the  shear  modulus 
has  very  little  influence  on  the  wrinkling  load,  f-or  the  case  of  a uniformly  loaded  beam 
it  is  found  that  the  wrinkling  load  varies  quadratically  with  the  pressure  and  linearly 
with  the  beam  cross-section  radius. 
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